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r(t) — e/e — l/\t\, and the system goes out of equilibrium at the time — i 
before the transition when r(—t) — t _1 (— t). After that time the state of the 
system essentially does not change until +i, when the rate r(t) becomes again 
equal to the relaxation rate, r(+t) — r~ 1 (+i). At i — T ^ , ^ 1+zu ^ fluctuations 

of the order parameter with wavelengths longer than | = (,(£) — Tq^ 1+zu ^ 
begin to grow exponentially while short wavelength fluctuations remain un- 
changed. The size of the correlated domains is given by | and the density of 

" 2 

vortices by where / ~ 10 in numerical experiments [7] as well as in 

analytic treatements of exactly soluble models [8] 

Theoretical and experimental studies of KZM have concentrated so far on 
thermal continuous phase transitions with the dynamics of the order param- 
eter governed by a phenomenological, irreversible time-dependent Ginzurg- 
Landau theory. Most attention has been devoted to the normal-superfluid tran- 
sition in A He [5], 3 He [6], superconductors [9, 10, 4], and, more recently, 
dilute Bose-Einstein condensates [11]. To date, truly microscopic approaches 
have been too complicated to extract useful predictions. 

A quench-induced quantum phase transition (QPT) at temperature T — 
must be treated in a microscopic way. It is a common wisdom that some prop- 
erties of a quantum transition can be obtained by an exact map from a ther- 
modynamic transition [12]: the correlation length in the ground state of the 
quantum system scales like £ ~ |e| _i/ and the gap between the ground state 
and the first excited state like dE ~ \e\~ zu . To study the dynamics of a quan- 
tum transition, which drives the system out of its ground state, we also need 
information about correlations in excited states, which is not provided by this 
map. As we will see below, the essence of the KZM, which is the competition 
between the transition rate and the timescale on which the system can react, 
is applicable to quantum transitions. However, the quantum scenario and, in 
particular, the interpretation of its results differ from the thermal case. The 
main reason is the reversibility of the quantum dynamics, as opposed to the 
dynamical irreversibility of the thermal critical dynamics. 

In this paper, we study the appearance of a nonzero current while the system 
is undergoing a (zero-temperature) quantum phase transition from the insulat- 
ing to the superfluid phase, providing details of a previously published study 
[13]. We also study the dynamics of the transition crossing the critical point 
in the opposite direction, namely from the superfluid to the insulator phase. 
When the transition is diabatic, this results in a finite dispersion in the number 
of atoms per site. We develop a microscopic, dynamical theory and suggest 
experiments to test our predictions. 
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2. Josephson junction arrays 

We consider a dilute gas of N atoms, trapped in a deep optical lattice. The 
lattice can be generated by counterpropagating lasers, which create a regular 
(possibly multidimensional) array of wells separated by potential barriers. The 
atoms are trapped on the bottom of each well, creating a lattice of weakly 
coupled Bose-Einstein condensates, the weak coupling being created by the 
tunneling of the atoms through the interwell barriers. When the height of the 
barrier is much larger than the "local" chemical potential of each condensate, 
the system can be described (in the "tight binding approximation") by the Bo- 
son Hubbard Model (BHM) Hamiltonian: 

h = -iE4^ + f £(«k) 2 , a) 

(k,i) k 

where g — An ah 2 /m is interaction strength, a is the interatomic s-wave scat- 
tering length and m the mass of the atoms. The term proportional to 7 in the 
BHM describes the atomic tunneling between the wells, while the term pro- 
portional to g(3 takes in account the interaction energy of the condensate atoms 
in each well. The coefficients (3 and 7 are integrals of the single particle bound 
state wave functions localized in each well [14]. 

We first consider the case when the total number of atoms N is commen- 
surate with the number of lattice sites, so that the average number of atoms 
per site n is an integer. The implications of a non-commensurate N will be 
discussed in Section 1.5.2. 

In the commensurate case the properties of the ground state of the BHM are 
governed by a dimensionless parameter 

G = =2 (2) 

It is possible to recognize four different regimes: 

■ n 2 <C G — the average interaction energy ngfi per atom can be ne- 
glected as compared to the hopping rate. The interaction term in the 
BHM (1) can be ignored and the ground state can be approximated by a 
quasi-coherent state 

\° s ) ~ (x>i) i°) ( 3 ) 

with the dispersion of the number of atoms in each site An ~ n 1 / 2 , 
and almost vanishing dispersion of the phase difference between nearest 
neighbor sites, 
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1 « G < n 2 — The on-site energy in the BHM is not anymore neg- 
ligible, but the atoms can still coherently tunnel from site to site. The 
interatomic interaction squeezes (i.e. reduces the fluctuations of) the 
number of atoms per site, An ~ G 1 / 4 . In this limit An is smaller than 
the coherent state fluctuations (~ n 1 / 2 ), yet is still greater than 1. Phases 
in different sites are correlated but their differences can have a significant 
dispersion which scales as 

A<j> ~ -!- ~ (T 1 / 4 . (5) 
An 

The phases become uncorrelated when G — > G+ ~ 0(1). At G = G c 
the system undergoes a quantum phase transition to an insulator (Mott) 
state characterized by a complete loss of phase coherence. 

G < G c — The ground state of the BHM can be very well approxi- 
mated by an incoherent superposition of localized states. At G = the 
ground state becomes 

\GS) = |n,n,n,...), (6) 



The nature of the transition can be captured by a simple argument. On the 
one hand, in the state (6), an atom in each site attempts to hop with a rate 7. 
The total hopping energy of the n atoms in each site is, therefore, nj. On the 
other hand, to move one atom to a nearest neighbor site costs an interaction 
energy ~ g(3. When G < G c — 0(1) the cost prevails over the gain provided 
by the hopping term, and the coherent tunneling is strongly suppressed. For 
any G < G c the dispersion of phase difference between nearest neighbor sites 
is maximal 

A<f> ~ 1 . (7) 

Such simple energetic arguments cannot prove that the critical point G = 
0(1) describes a phase transition and not a crossover. We note, therefore, that 
the zero temperature model (1) in d dimensions can be exactly mapped on the 
classical X-Y model at finite temperature in d+ 1 dimensions [12]. The X — Y 
model has a continuous thermal phase transition. Thermal correlators in the X- 
Y model can be mapped back to quantum correlators in our quantum system. 
We will use this mapping in Section 1.6. 
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3. The quantum phase model 

For G<n 2 and n > 1 the boson Hubbard model (1) can be approximated 
by a quantum phase model (QPM) [15]: 

H = E 9 4tf ~ 7 E cos (4 - k) (8) 
k (k,l) 

with the last sum running over the nearest-neighbor sites. The phase <f>i and 
the number of atoms ni in each site of the junction are (with some important 
caveats [16]) non-commuting conjugate observables [ni,4>i] — i, and, in the 
^-representation, 



ni = n- i— , (9) 
o<Pl 

4>i = 4>i- (10) 

Therefore, cf>i and ni play the role of coordinate and conjugate momentum, and 
satisfy the Heisenberg uncertainty relation. The insulator phase is character- 
ized by large quantum phase fluctuations in the ground state, which destroy the 
long-range order among sites. The ID Hamiltonian Eq.(8) exhibits a continu- 
ous Mott phase transition at the critical value G c — 0.617 [12]. 

The QPM becomes equivalent to the BHM (1) in the limit of large number 
of atoms per site, n >• 1, and for G <n 2 [17] (when G 3> n 2 the Schredinger 
equation contains extra cos 2(^ — fy) terms which have been neglected here). 
These conditions are well satisfied in current experiments where ~ 10 5 con- 
densate atoms are trapped in ~ 10 3 wells created by a one dimensional optical 
lattice. Therefore, the QPM provides an unified model to study a nonequilib- 
rium QPT in different systems, and, at the same time, allows the development 
of simple approximate methods and of a more transparent physical intuition of 
the process. 

In what follows we study the microscopic dynamics of quantum phase tran- 
sitions in a one-dimensional chain with periodic boundary conditions. The 
dynamics is governed by the QPM Eq.(8) which after rescaling time 

t -+ gfit (11) 
can be written as a dimensionless Schredinger equation 

4*=-5?M*" G ^ cosWt "* ,) *- <12) 

k t k (k,l) 
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The ground state below the transition (G < G c ) is close to the localized state 
(6) with all sites occupied by the same number of atoms n. In the phase rep- 
resentation this state is described by a uniform wavefunction, ^(4>i) — const, 
where all phase differences between nearest neighbor sites have the maximal 
dispersion of A<p ~ 1. Above the transition (G > G c ) the ground state 
is a number squeezed state which continuously tends to a Fock state when 
G — > G+, and to a coherent state for G 3> n 2 . Therefore, when G 3> G c and 
G«n 2 one can describe the low energy part of the spectrum of Eq.(12) in a 
harmonic approximation [18] 

# = + f £(^) 2 *. d3) 

m 1 k a< Pk 1 <fcji> 

This equation is diagonalized by normal modes numbered by a lattice momen- 
tum 

fiE{-N s + l,...,+N s }, (14) 

* — Ylfi ^ni^n)- H ere N s is the number of lattice sites. There is one zero 
mode $o ~ X)i All other modes $ M have nonzero frequencies y^G, 



.d 

'at*" 



2 m. * 2 " " 



1 — cos 



7T/i 



7/i = 2 

We remark that all frequencies scale as 



(15) 
(16) 

(17) 



and phase dispersions in the ground states of the harmonic oscillators (15) scale 
as 

A% = ( 7i uG)- 1/4 ~ G- 1 / 4 . (18) 
We will use these scalings in the following. 

4. Linear Quench 

To quench the system trough the transition point in a timescale tq, we lin- 
early ramp the control parameter G in Eq.(12) as 



G(t) 



r Q 



(19) 
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The time t runs from to tmax = TQG max with G max > G c . 

The initial state is the localized ground state at G — 0, see Eq.(6). The goal 
is to preserve memory of the phase disorder characteristic of the low G Mott 
insulator phase into the large G phase ordered superfluid phase. We begin with 
two technically trivial limits, tq — > oo and tq — > 0. 

4.1 Adiabatic transition: tq — >■ oo 

In the limit tq — > oo the transition is adiabatic and the state of the system 
follows its instantaneous ground state from the initial localized state |n, . . . , n) 
at G — to a phase squeezed state with all the nearest neighbor phase differ- 
ences close to zero, A<f> « 0, at G max 3> G c . The ground state at G m3x 3> G c 
is well described by the harmonic oscillators (15). Therefore, the final disper- 
sion of phase differences in the adiabatic transition is 



The phase dispersion in the final squeezed state can be probed with interference 
measurements. 

4.2 Instantaneous transition: r Q — > 

In the opposite limit of an instantaneous transition, when tq — > 0, the sys- 
tem is not able to adjust its quantum state to the changing G(t), and it remains 
in the initial localized state (6) till G max . The final dispersion of the phase 
differences in the uniform Fock state, *(</>/) = const, is 



The relative phases are random and we expect no interference fringes. 

4.3 Diabatic transition 

In the following we consider diabatic transitions in intermediate regimes. In 
Fig.(l) we illustratively show inverse gap r = 1/dE of the system as a func- 
tion of the control parameter G(t) . The inverse gap is a time scale on which the 
quantum system can react to changes of the external parameter G(t). It plays 
the same role as the relaxation time in the thermal KZM. The inverse gap van- 
ishes at the transition point G — G c . For a generic quench the evolution of the 
system is approximately adiabatic when G is far from G c , and impulse when G 
is close to G c ; compare with Fig.(2) where we consider a 3-site periodic lattice. 
This also happens in a thermal phase transition [2]. As G(t) increases starting 
from 0, the system initially remains in its instantaneous ground state, which is 
an incoherent superposition of localized states, until the point G_ < G c where 
the transition ceases to be adiabatic. Between G_ and a certain G > G c the 




(20) 



(21) 
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Figure 1. The inverse gap of the system as a function of the control parameter G{t). The 
transition point is at G c — t c / tq . The transition is impulse when G- < G < G. The time 
runs from to t ma x — TQGmax ■ 

evolution is impulse: to a first approximation the state of the system does not 
change. At the point G(t) — G the system is still described by an incoherent 
superposition of localized states with a large phase dispersion A(f> ~ 1. When 
G(t) > G the transition becomes adiabatic again. 
We consider two diabatic regimes 

■ the gaussian regime, which is reached when G c <C G. In this regime we 
can use the gaussian approximation (15) when G(t) ~ G. 

■ the critical regime, reached when G ~ G c . In this case the gaussian ap- 
proximation (15) breaks down, and we need to use exact renormalization 
group critical exponents. 



5. The gaussian regime 

With the linear quench (19) the periods of the oscillators decay with time 

as: 

(my*- M 1 ". <22) 



\GuiaxJ \ t J 

As shown in Appendix, the instantaneous time scale of the transition is 
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£=t. (23) 

dt 

This time scale becomes comparable to the periods of the short wavelength 
oscillators with « N s at the time t when G -1 / 2 ~ G/G, see Fig.(2). The 
solution of this equation gives: 

t ~ r^ 3 , (24) 
G = G(t) ~ Tq 2/3 . (25) 

For i > t the short wavelength oscillators respond adiabatically to the quench. 
Oscillators with longer wavelengths become adiabatic later, but still at a time 
scaling as t ~ Tq 3 . The system is in the gaussian regime, G c <C G, when 

tq < 1 . (26) 

The long wavelength oscillators become adiabatic at G^^i ~ {N s /tq)~ 2 ^. 
To remain in the gaussian regime, we want G|^| = i <C n 2 . This condition trans- 
lates into 

N s < tq n 3 . (27) 

For, say, n ~ 10 3 atoms per site this condition is satisfied in a wide range of 
lattice sizes N s and quench times tq. 

We now estimate the final dispersion A(f> of phase differences between the 
nearest neighbor lattice sites. At G the dispersion is A<f> ~ 1 like in the initial 
Mott insulator state (6). For G > G the evolution of the harmonic oscillators 
is adiabatic. There is no mixing between the eigenstates of any oscillator fj, 
(15). The dispersion of the phase in any given eigenstate of the oscillator 
(15) scales as G -1 / 4 . In the adiabatic evolution after G the dispersion of ^ 
shrinks like [G/G^)] 1 / 4 . A<p behave in the same way. The final dispersion at 

Grriax iS 



- r^G-V 4 . (28) 

"max J 

In a ID periodic lattice this A<f> translates into a dispersion AL of the "angular 
momentum" 

L - i ^2(ajai +1 - a\ +1 a{) , (29) 
i 

which depends on the nearest neighbor A(j>, through the formula 
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Figure 2. The inverse of the gap T — 1 /dE as a function of G for a 3-site periodic Hubbard 
model with n — 6 atoms per site (dots). The solid line is the best G -1 / 2 fit. The dashed lines 
show the transition time G/G — t — GrQ/Gma* for G ma x = 10 and tq — 4.5, 2.1, 0.9 (from 
left to right). The crossings between 1/dE and the dashed lines define G for the different tq. 



AL ~ VN s nA(j) (30) 

valid for a small A<f>. This AL may translate into a dispersion of the wind- 
ing number y/Ng A<f> when the atoms after the quench are somehow forced to 
condense. 

Figures 3 and 4 show results of numerical simulations of the boson Hubbard 
model with a 3-site periodic lattice with a total of 300 atoms. In Fig.3 we show 
that A<p remains the same as in the initial Fock state for G < G T q<i • hi Figure 

1/3 

4 we verify the scaling t ~ Tq . 

5.1 Soft modes 

In the argument of this Section we have neglected so far that different nor- 
mal modes have different frequencies, see Eq.(16). We used only the general 
property that all frequencies, except for the zero mode, scale with G as G 1//2 . 
However, it is clear that long wavelength soft phonon modes will become adi- 
abatic later than short wavelength modes. We can even define a G^ for each 
mode fx as a G when the frequency of the mode becomes comparable to the 
transition rate, 
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Figure 3. The normalized dispersion dl — AL/AL TQ _>o as a function of 5 = G(t) /G max . 
We show results of exact numerical simulations of the 3-site Hubbard model for 7 different 
tq - 0.016, 0.008, 0.004, 0.002, 0.001, 0.0005, 0.00025 (from left to right) with G ma x = 
1000 and n — 100 atoms per site, dl deviates from 1 (and hits the bottom of the figure) at 
G ~ Tq 2 ^ 3 : compare Fig.(4). 



\fh£^S- (3D 
Solving this equation with respect to G we obtain 

* % 1/3 t q 2/3 . (32) 

After Gn the dispersion of the phase ^ shrinks like G - 1 / 4 and at a G max it 
becomes 



A*„ c c ^\ 1/6 G^ , (33) 

compare with Eq.(28). The dependence on 7^ is much softened by the —1/12 
exponent. 7^ ranges from around 7r 2 /iVf for fj, — 1 to 2 for fj, — N s /2, 

—1/12 

compare Eq.(16). For, say, N s — 100 lattice sites in ID the prefactor 7^ ' 
changes by less than a factor of 2. This justifies a posteriori the implicit ap- 
proximation in our argument that 7^ ~ 0(1). 
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Figure 4. The normalized dispersion = AL/AL TQ _>o as a function of s — G(t) /G with 
G ~ Tq 2 ^ 3 . This figure displays all the seven plots from Fig.3. The plots in Fig. 3 are in the top 
left corner of this figure. The fact that the seven plots sit on top of each other proves that the 
evolution of AL depends on time t through the combination t/i with i ~ Tq 3 . 

5.2 Commensurate versus non-commensurate 

So far we have assumed that the number of atoms N is commensurate 
with the number of sites. In a rigorous mathematical sense, when N is non- 
commensurate, there is no phase transition. The system is always in a super- 
fluid phase. It is easy to see why. Suppose that we have n atoms at each site 
and we add 1 atom to the system in order to make N non-commensurate. The 
extra atom can hop between sites with a rate 7, since there is no interaction 
energy cost to suppress this hopping. The ground state for G -C 1 is 



0.8 



0.6 — 



\GS) ~ \E a k) l«> «>«>■■■>■ (34) 

The bulk of commensurate atoms are in an insulating state but the extra non- 
commensurate atom is superfluid even for this very small G. 

In Fig.5 we show results of the same numerical simulation as in Fig.3 but 
with N — 300 + 1 atoms or n — 100 + ^ 3> 1 atoms per site. The plots in 
Fig.5 are essentially identical with the plots in Fig.3. This demonstrates that, at 
least in the gaussian regime, the transition is dominated by the commensurate 
bulk. 
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Figure 5. The normalized dispersion dl — AL/AL TQ _>o as a function of 5 = G(t) /G max . 
We show results of exact numerical simulations of the 3-site Hubbard model for 7 different 
tq - 0.016, 0.008, 0.004, 0.002, 0.001, 0.0005, 0.00025 (from left to right) with G max = 
1000 and n — 100 + | atoms per site. These plots are indistinguishable from plots in Fig.3 
obtained for an integer n — 100. 



6. The critical regime 

So far we have considered the case tq <C 1 with G well inside the harmonic 
regime. In slower transitions, when G ~ G c , the harmonic approximation 
cannot be applied. The slow transitions (tq <C 1) probe the critical behavior 
of the system close to G c . 

■ In the periodic 3-site array there is no phase transition but, rather, a 
crossover: the gap dE is minimal at the crossover point G c ~ 0.5 but 
it does not vanish there, compare Figs.(l,2). For tq ^> 1 the transi- 
tion is adiabatic and the system follows its instantaneous ground state. 
The phase difference dispersion at G max >> G c is the dispersion in the 
ground state of the harmonic oscillators (20). A thermal crossover tran- 
sition was studied in [21]. 

■ In a large 1 D array with a commensurate density of atoms there is a 
phase transition and not a crossover [12, 19]. The gap dE vanishes 
at G — G c . The quantum zero temperature transition translates into 
Berezinskii-Kosterlitz-Thouless transition which does not have a well 
defined correlation length in the ordered phase - correlations decay alge- 



braically with a distance on a lattice. The scaling argument alone is not 
enough to make predictions of A<f> in this case. 

In a large 2 D array with a commensurate density of atoms there is a 
second order phase transition. Close to the critical point, G « G c , the 
energy gap scales as 



dE ~ (G - G c y , (35) 

where v = 2/3 is an exact renormalization group critical exponent of 
the X-Y model in 2 + 1 D. 1 With a linearized 



G(t) - G c = — (36) 
T Q 

the gap AE ~ (G — G c ) 2 / 3 = (*/tq) 2 / 3 becomes equal to the quench 
rate G/(G-G C ) = 1/t at 

G^Gc + ^§. (37) 
T Q 

Before G the system is in an incoherent superposition of localized states. 
After G the evolution is adiabatic, and with increasing G the phase dis- 
persion shrinks together with the phase widths of the system eigenstates. 
The dispersion measured at G max ^> G c , where we can use the harmonic 
approximation (13,18), is 



"max umax 




This formula is consistent with Eq.(20) because G c — 0(1). Here we 
again use the fact that the phase dispersion of harmonic oscillator eigen- 
states shrinks like G -1 / 4 . The critical behavior is more transparent in a 
phase dispersion offset with a phase dispersion for large tq, 



- A^oo ~ 3/5 . (39) 

t^max Tq 

The critical exponent u — 2/3 determines the exponent of 3/5. 



In a large 3 D array there is a second order phase transition with a gaus- 
sian critical exponent of v — 1/2 characteristic for the X-Y model in 



Non Equilibrium Mott Transition 



15 



3 + 1 D. In 3 D the results of the gaussian regime are valid also in the 
critical regime. The Mott transition has been recently observed in 1 D 
[22], and in 3 D [23] . 

7. From superfluid to insulating phase 

So far we have focused on transitions from the insulating to the superfluid 
phase where the non-equilibrium dynamics of the transition leaves imprint on 
the dispersion of phase difference between sites A(j>. The phase dispersion can 
be measured by interference techniques [20]. In this Section we will reverse 
the direction of the transition to go from the superfluid to the insulating phase. 
We will use the dispersion in the number of atoms per site An as a stable record 
of the nonequilibrium dynamics of the transition. In the insulating phase the 
number of atoms in each site is conserved. 

Again we linearly ramp G from a G max ( 1 <C G max <C n 2 ) at t — 
— T^Gmax down to G — at t = 0, 

G(t) = — . (40) 
We can distinguish several regimes for different values of tq. 

■ Instantaneous transition with tq — > 0. The system remains in the 
initial ground state at G m ax which is a squeezed quasi-coherent state 
with a dispersion 

An TQ ^ ^ GUL • (41) 

■ Adiabatic limit of tq — > oo. The state of the system follows the in- 
stantaneous ground state. The system ends in the Mott insulator state (6) 
with a dispersion 

An TQ ^oo = . (42) 

■ Diabatic transition with tq <C 1 - the gaussian regime. The char- 
acteristic frequency in the harmonic regime scales like G 1 / 2 , compare 
Eq.(17). This frequency becomes equal to the transition rate G/G at a i 

—2/3 

when G ~ Tq . Before t the system follows its ground state which is 

a squeezed state with An ~ G 1 / 4 . At i when An ~ G 1 / 4 the evolution 
becomes impulse, the state does not change with the decreasing G(t) 
any more. The system arrives at G = with 



An TQ<<1 ~ G 1 / 4 ~ tq 



(43) 
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0.0 1.0 2.0 

Figure 6. Logarithm of the dispersion of the number of atoms per site Y — log 10 An as a 
function of X — log 10 tq . The dots are averages over many numerical simulations of the 3-site 
periodic Hubbard model (1) for n — 100 atoms per site. The solid line is the best linear fit with 
a slope of —0.164 ± 0.002. This slope compares well with the prediction of —1/6 — —0.167 
inEq.(43). 

This scaling applies when G <C G max or tq Gm^ 2 . In the opposite 
case the transition is instantaneous. 

In Fig.6 we show An as a function of tq in a 3-site periodic Hubbard 
model with n — 100 atoms per site. The best fit to the numerical re- 
sults gives an exponent of —0.164 ± 0.002 which is consistent with the 
prediction of -1/6 = 0.167 in Eq.(43). 

■ Diabatic transition with tq 1 - the critical regime. For a 2D or 3D 

array the gap dE depends on a distance from the critical point \G C \ as 
dE ~ \G — G c \ v with v — | in 2D and v — \ in 3D. The dispersion An 
in the ground state at G c is An c 0.5 and close to G c the dispersion can 
be linearized as An— An c ~ G—G c — —t/TQ. The evolution is impulse 
between two times, —i and +t, when the transition rate G/(G — G c ) — 
l/\t\ equals the gap dE ~ \G — G^. This happens when \G — G c \ ~ 
T Q 1 ^ 1+ ^ ■ ^ ne system adiabatically follows its instantaneous ground 
state until — i when the dispersion in the ground state is An_ — An c ~ 
+\G — G c \. Between — i and +t the evolution is impulse, in the first 
approximation the state of the system does not change. The dispersion 
in the instantaneous ground state at +i is An + — An c ~ — \G — G c \ but 
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the state of the system still has the dispersion An_ > An+. This extra 
dispersion above that in the instantaneous ground state is 

An- - An + ~ r^ 1 ^ 7 . (44) 

This is dispersion at +i but we want to know the final dispersion when 
G = 0. 

In the insulating phase the An of the ground state strongly depends on 
G, it ranges from An = at G — to An c « 0.5 at G — G c « 1. In 
contrast, An of excited states remains nonzero even at G — 4 " because 
an arbitrarily small but nonzero tunneling rate G removes degeneracy of 
states with different occupation numbers. All excited states at G — + 
have nonzero An. Their An depends on G but, in first approximation, 
this dependence can be neglected as compared to that in the ground state. 
In short, after +i the An of the ground state shrinks down to while the 
An's of the excited states does not change. The state of the system at 
t — +i is a superposition of the ground state and the excited states. The 
adiabatic evolution after +t ends in a state with a An that is a linear 
function of the difference in Eq.(44), 

An TQ »i ~ . (45) 

This final An depends both on the transition time tq and on the critical 
exponent v. 

8. Conclusion 

After crossing the transition, the system does not have a definite angular 
momentum or definite phase differences between lattice sites, but it is in a 
quantum superposition of states with different Acf>. Either a measurement or 
decoherence [24] are needed to convert this coherent superposition into a mix- 
ture of states, each with definite current. This is an important difference with 
respect to the case of a thermal KZM, where A<f> would describe the disper- 
sion in an ensemble of different possible classical outcomes. Moreover, the 
quantum dynamics is reversible, while in the thermal case the characteristic 
lenghtscale £ is frozen after the symmetry breaking transition is completed. £ 
is a record of the transition rate tq : manipulations with e do not change the 
winding number as long as the system remains in the symmetry broken phase 
with e > 0. In contrast, even after the diabatic quantum transition is com- 
pleted one can change A<p TQ . Adiabatic variations of G away from G max are 
accompanied by changes in A<f), 



A^A-MGnWG) 1 / 4 



(46) 
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In conclusion, we have predicted the phase dispersion after a diabatic insulator- 
superfluid quantum phase transition in an array of weakly coupled Bose-Einstein 
condensates or Josephson junctions. This theory is a quantum counterpart of 
the Kibble-Zurek mechanism for topological defect formation in classical ther- 
mal phase transitions. Our predictions can be tested experimentally with a 
Bose-Einstein condensate trapped in a periodic potential [20, 23]. In this case 
our Eq.(28) predicts a dispersion of the phases A<f> ~ G^J^Tq 1//6 for a transi- 
tion to the superfiuid phase and a dispersion of the number of atoms per lattice 

site of An ~ Tq 1 ^ 6 for a transition to the insulating phase. The dispersions can 
be directly measured by, respectively, interference and phase contrast imaging 
techniques. 
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Appendix 

Here we justify the choice of ^ /G, and not, say, ^ — Tq 1 , as the transition rate. The 
wavefunction in Eq.(15) can be expanded in eigenstates of the harmonic oscillator, 

oo 

where U m [x] are eigenstates of a dimensionless harmonic oscillator 

+ r Um[x] = ( m + D Um[x] ■ (A - 2) 

Substitution of the expansion (A.l) to the Schredinger equation (15) gives 

oo 

El 
- 

oo 

E> 

m=0 

Herea; = $ M /(7 M G) 1/4 . We can eliminate the first term on the right hand side by a redefinition 
of the amplitude 

(*) — B m (t) exp 
This is just a change of phase factors. After projection on U n [x] we get 



(7mG) 1/4 



(A.1) 




dt - 4\G J 
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f; J? m e i(n - ro) /o j dxxU n ^- . (A.4) 

771=0 ^ 

(dG \ 
-^f- J sets the rate 

at which the B n is forced to change. However, the right hand side oscillates with frequencies 
which are multiplicities of the instantaneous frequency of the mode \/"/ij,G(t). The effect of 
the right hand side on B n averages to zero when 

/ dG \ 

i^-f) « (A.5) 
With this condition different amplitudes B n effectively decouple and the evolution is adiabatic. 
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